Variability of the innate immune response is globally constrained by transcriptional bursting

Transcription of almost all mammalian genes occurs in stochastic bursts, however the fundamental control mechanisms that allow appropriate single-cell responses remain unresolved. Here we utilise single cell genomics data and stochastic models of transcription to perform global analysis of the toll-like receptor (TLR)-induced gene expression variability. Based on analysis of more than 2000 TLR-response genes across multiple experimental conditions we demonstrate that the single-cell, gene-by-gene expression variability can be empirically described by a linear function of the population mean. We show that response heterogeneity of individual genes can be characterised by the slope of the mean-variance line, which captures how cells respond to stimulus and provides insight into evolutionary differences between species. We further demonstrate that linear relationships theoretically determine the underlying transcriptional bursting kinetics, revealing different regulatory modes of TLR response heterogeneity. Stochastic modelling of temporal scRNA-seq count distributions demonstrates that increased response variability is associated with larger and more frequent transcriptional bursts, which emerge via increased complexity of transcriptional regulatory networks between genes and different species. Overall, we provide a methodology relying on inference of empirical mean-variance relationships from single cell data and new insights into control of innate immune response variability

Heat shock response regulates stimulus-specificity and sensitivity of the pro-inflammatory NF-κB signalling

From Springer Nature via Jisc Publications RouterHistory: received 2020-01-24, accepted 2020-04-16, registration 2020-04-16, pub-electronic 2020-05-24, online 2020-05-24, collection 2020-12Publication status: PublishedFunder: Narodowe Centrum Nauki; doi: http://dx.doi.org/10.13039/501100004281; Grant(s): 2016/23/B/ST6/03455, 2015/19/B/ST7/02984, 2016/21/B/ST7/02241Funder: Politechnika Śląska; doi: http://dx.doi.org/10.13039/501100007835; Grant(s): 02/010/BK19/0143Funder: Biotechnology and Biological Sciences Research Council; doi: http://dx.doi.org/10.13039/501100000268; Grant(s): BB/K003097/1Funder: Politechnika Poznańska; doi: http://dx.doi.org/10.13039/501100004239; Grant(s): 825/RN2/RR4/2018Abstract: Background: Ability to adapt to temperature changes trough the Heat Shock Response (HSR) pathways is one of the most fundamental and clinically relevant cellular response systems. Heat Shock (HS) affects the signalling and gene expression responses of the Nuclear Factor κB (NF-κB) transcription factor, a critical regulator of proliferation and inflammation, however, our quantitative understanding of how cells sense and adapt to temperature changes is limited. Methods: We used live-cell time-lapse microscopy and mathematical modelling to understand the signalling of the NF-κB system in the human MCF7 breast adenocarcinoma cells in response to pro-inflammatory Interleukin 1β (IL1β) and Tumour Necrosis Factor α (TNFα) cytokines, following exposure to a 37–43 °C range of physiological and clinical temperatures. Results: We show that exposure to 43 °C 1 h HS inhibits the immediate NF-κB signalling response to TNFα and IL1β stimulation although uptake of cytokines is not impaired. Within 4 h after HS treatment IL1β-induced NF-κB responses return to normal levels, but the recovery of the TNFα-induced responses is still affected. Using siRNA knock-down of Heat Shock Factor 1 (HSF1) we show that this stimulus-specificity is conferred via the Inhibitory κB kinase (IKK) signalosome where HSF1-dependent feedback regulates TNFα, but not IL1β-mediated IKK recovery post HS. Furthermore, we demonstrate that through the temperature-dependent denaturation and recovery of IKK, TNFα and IL1β-mediated signalling exhibit different temperature sensitivity and adaptation to repeated HS when exposed to a 37–43 °C temperature range. Specifically, IL1β-mediated NF-κB responses are more robust to temperature changes in comparison to those induced by TNFα treatment. Conclusions: We demonstrate that the kinetics of the NF-κB system following temperature stress is cytokine specific and exhibit differential adaptation to temperature changes. We propose that this differential temperature sensitivity is mediated via the IKK signalosome, which acts as a bona fide temperature sensor trough the HSR cross-talk. This novel quantitative understanding of NF-κB and HSR interactions is fundamentally important for the potential optimization of therapeutic hyperthermia protocols. D-ESB-DSZSgbr1c4s2oumkVideo Abstrac

Live-cell imaging reveals single-cell and population-level infection strategies of Listeria monocytogenes in macrophages

Pathogens have developed intricate strategies to overcome the host’s innate immune responses. In this paper we use live-cell microscopy with a single bacterium resolution to follow in real time interactions between the food-borne pathogen L. monocytogenes and host macrophages, a key event controlling the infection in vivo. We demonstrate that infection results in heterogeneous outcomes, with only a subset of bacteria able to establish a replicative invasion of macrophages. The fate of individual bacteria in the same host cell was independent from the host cell and non-cooperative, being independent from co-infecting bacteria. A higher multiplicity of infection resulted in a reduced probability of replication of the overall bacterial population. By use of internalisation assays and conditional probabilities to mathematically describe the two-stage invasion process, we demonstrate that the higher MOI compromises the ability of macrophages to phagocytose bacteria. We found that the rate of phagocytosis is mediated via the secreted Listeriolysin toxin (LLO), while the probability of replication of intracellular bacteria remained constant. Using strains expressing fluorescent reporters to follow transcription of either the LLO-encoding hly or actA genes, we show that replicative bacteria exhibited higher PrfA regulon expression in comparison to those bacteria that did not replicate, however elevated PrfA expression per se was not sufficient to increase the probability of replication. Overall, this demonstrates a new role for the population-level, but not single cell, PrfA-mediated activity to regulate outcomes of host pathogen interactions

Integration of kinase and calcium signaling at the level of chromatin underlies inducible gene activation in T cells

TCR signaling pathways cooperate to activate the inducible transcription factors NF-κB, NFAT, and AP-1. In this study, using the calcium ionophore ionomycin and/or PMA on Jurkat T cells, we show that the gene expression program associated with activation of TCR signaling is closely related to specific chromatin landscapes. We find that calcium and kinase signaling cooperate to induce chromatin remodeling at ∼2100 chromatin regions, which demonstrate enriched binding motifs for inducible factors and correlate with target gene expression. We found that these regions typically function as inducible enhancers. Many of these elements contain composite NFAT/AP-1 sites, which typically support cooperative binding, thus further reinforcing the need for cooperation between calcium and kinase signaling in the activation of genes in T cells. In contrast, treatment with PMA or ionomycin alone induces chromatin remodeling at far fewer regions (∼600 and ∼350, respectively), which mostly represent a subset of those induced by costimulation. This suggests that the integration of TCR signaling largely occurs at the level of chromatin, which we propose plays a crucial role in regulating T cell activation

Interactions among oscillatory pathways in NF-kappa B signaling

<p>Abstract</p> <p>Background</p> <p>Sustained stimulation with tumour necrosis factor alpha (TNF-alpha) induces substantial oscillations—observed at both the single cell and population levels—in the nuclear factor kappa B (NF-kappa B) system. Although the mechanism has not yet been elucidated fully, a core system has been identified consisting of a negative feedback loop involving NF-kappa B (RelA:p50 hetero-dimer) and its inhibitor I-kappa B-alpha. Many authors have suggested that this core oscillator should couple to other oscillatory pathways.</p> <p>Results</p> <p>First we analyse single-cell data from experiments in which the NF-kappa B system is forced by short trains of strong pulses of TNF-alpha. Power spectra of the ratio of nuclear-to-cytoplasmic concentration of NF-kappa B suggest that the cells' responses are entrained by the pulsing frequency. Using a recent model of the NF-kappa B system due to Caroline Horton, we carried out extensive numerical simulations to analyze the response frequencies induced by trains of pulses of TNF-alpha stimulation having a wide range of frequencies and amplitudes. These studies suggest that for sufficiently weak stimulation, various nonlinear resonances should be observable. To explore further the possibility of probing alternative feedback mechanisms, we also coupled the model to sinusoidal signals with a wide range of strengths and frequencies. Our results show that, at least in simulation, frequencies other than those of the forcing and the main NF-kappa B oscillator can be excited via sub- and superharmonic resonance, producing quasiperiodic and even chaotic dynamics.</p> <p>Conclusions</p> <p>Our numerical results suggest that the entrainment phenomena observed in pulse-stimulated experiments is a consequence of the high intensity of the stimulation. Computational studies based on current models suggest that resonant interactions between periodic pulsatile forcing and the system's natural frequencies may become evident for sufficiently weak stimulation. Further simulations suggest that the nonlinearities of the NF-kappa B feedback oscillator mean that even sinusoidally modulated forcing can induce a rich variety of nonlinear interactions.</p

Modeling stochasticity in gene regulation

Intrinsic stochasticity plays an essential role in gene regulation because of the small number of involved molecules of DNA, mRNA and protein of a given species. To better understand this phenomenon, small gene regulatory systems are mathematically modeled as systems of coupled chemical reactions, but the existing exact description utilizing a Chapman-Kolmogorov equation or simulation algorithms is limited and inefficient. The present work introduces a much more efficient yet accurate modeling approach, which allows analyzing stochasticity in the system in terms of the underlying distribution function. The novel modeling approach is motivated by the analysis of a single gene regulatory module with three sources of stochasticity: intermittent gene activity, mRNA transcription/decay and protein translation/decay noise. Although the corresponding Chapman Kolmogorov equation cannot be solved when a large number of molecules are considered, it is used to analytically derive the first two moments of the underlying distribution function. The mRNA and protein variance is found decomposable into additive terms resulting from the respective sources of stochasticity, which allow quantifying their significance in the process. The variance decomposition is asserted by constructing two approximations that establish a novel modeling approach: First, the continuous approximation, which considers only the stochasticity due to the intermittent gene activity. Second, the mixed approximation, which in addition attributes stochasticity to the mRNA transcription/decay process. Introduced approximations yield systems of first order partial differential equations for the underlying distribution function, which can be efficiently solved using developed numerical methods. Single cell simulations and numerical two-dimensional mRNA-protein stationary distribution functions are presented to confirm accuracy of introduced models. Further simplifications in the model allow considering regulation of the two- (possibly three-) gene systems for which two-dimensional protein-protein distributions are calculated. Finally, the assumption that gene activity is due to the binding and dissociation of a single regulatory molecule is relaxed. Based on the gene expression data, the models developed are applied to hypothesize the existence of a sequential activation mechanism of NF-kappaB dependent genes important in cell survival and inflammation. Future applications include analysis of small genetic networks, which are being currently engineered based on the prokaryotic and eukaryotic components

Application of Sensitivity Analysis to Discover Potential Molecular Drug Targets

Mathematical modeling of signaling pathways and regulatory networks has been supporting experimental research for some time now. Sensitivity analysis, aimed at finding model parameters whose changes yield significantly altered cellular responses, is an important part of modeling work. However, sensitivity methods are often directly transplanted from analysis of technical systems, and thus, they may not serve the purposes of analysis of biological systems. This paper presents a novel sensitivity analysis method that is particularly suited to the task of searching for potential molecular drug targets in signaling pathways. Using two sample models of pathways, p53/Mdm2 regulatory module and IFN- [Formula: see text]-induced JAK/STAT signaling pathway, we show that the method leads to biologically relevant conclusions, identifying processes suitable for targeted pharmacological inhibition, represented by the reduction of kinetic parameter values. That, in turn, facilitates subsequent search for active drug components